Contemporary Glacier Processes and Global Change: Recent Observations from Kaskawulsh Glacier and the Donjek Range, St. Elias Mountains

With an extensive ice cover and rich display of glacier behaviour, the St. Elias Mountains continue to be an enviable natural laboratory for glaciological research. Recent work has been motivated in part by the magnitude and pace of observed glacier change in this area, which is so ice-rich that ice loss has a measurable impact on global sea level. Both detection and attribution of these changes, as well as investigations into fundamental glacier processes, have been central themes in projects initiated within the last decade and based at the Kluane Lake Research Station. The scientific objectives of these projects are (1) to quantify recent area and volume changes of Kaskawulsh Glacier and place them in historical perspective, (2) to investigate the regional variability of glacier response to climate and the modulating influence of ice dynamics, and (3) to characterize the hydromechanical controls on glacier sliding. A wide range of methods is being used, from ground-based manual measurements to space-based remote sensing. The observations to date show glaciers out of equilibrium, with significant ongoing changes to glacier area, volume, and dynamics. Computer models are being used to generalize these results, and to identify the processes most critical to our understanding of the coupled glacier-climate system.Grâce à leur importante couverture de glace et au riche étalage de comportement des glaciers, les monts St. Elias continuent de servir de laboratoire naturel enviable pour la recherche glaciologique. Des études récentes ont été motivées, en partie, par la magnitude et la vitesse des changements observés dans les glaciers de l’endroit, qui sont riches en glace au point que la perte de glace a une incidence mesurable sur le niveau général de la mer. La détection et l’attribution de ces changements de même que les recherches à l’égard des processus des glaciers ont servi de thème central à des projets qui ont été mis en oeuvre au cours de la dernière décennie à la station de recherche du lac Kluane. Les objectifs scientifiques de ces projets consistent (1) à quantifier les changements récents relativement à l’aire et au volume du glacier Kaskawulsh, puis à les mettre dans une perspective historique, (2) à faire enquête sur la variabilité générale de la réaction du glacier vis-à-vis du climat et de l’influence modulatrice de la dynamique de la glace, et (3) à caractériser le contrôle hydromécanique par rapport au glissement du glacier. Une vaste gamme de méthodes est employée pour parvenir à ces fins, allant des mesures manuelles sur le terrain à la télédétection spatiale. Jusqu’à maintenant, les observations indiquent que les glaciers ne sont pas en équilibre et que d’importants changements se produisent quant à l’aire, au volume et à la dynamique du glacier. Des modèles informatiques sont utilisés pour généraliser ces résultats ainsi que pour cerner les processus les plus critiques à notre compréhension du système couplé glacier-climat

Modeling Sediment Transport in Ice-Walled Subglacial Channels and Its Implications for Esker Formation and Proglacial Sediment Yields

Sediment yields from glacierized basins are used to quantify erosion rates on seasonal to decadal timescales as well as conditions at the glacier bed, and eskers hold valuable information about past subglacial hydraulic conditions in their spatial organization, geometry, and sedimentary structures. Ultimately, eskers are a record of past glacio‐fluvial sediment transport, but there is currently no physical model for this process. We develop a 1‐D model of morphodynamics in semicircular bedrock‐floored subglacial channels. We adapt a sediment conservation law developed for mixed alluvial‐bedrock conditions to subglacial channels. Channel evolution is a function of melt opening by viscous heat dissipation from flowing water and creep closure of the overlying ice, to which we add the closure or enlargement due to sediment deposition or removal, respectively. We apply the model to an idealized land‐terminating glacier and find that temporary sediment accumulation in the vicinity of the terminus, or the formation of an incipient esker, is inherent to the dynamics of the channelized water flow. The alluviation of the bed combined with the pressurized channel flow produces unexpected patterns of sediment evacuation: We show that the direction of hysteresis between sediment and water discharge is not necessarily linked to a supply‐ or transport‐limited system, as has been hypothesized for proglacial sediment yields. We also find that the deposition of an incipient esker is a function of a compromise between water discharge and sediment supply, but perhaps more importantly, ice‐surface slope and the temporal pattern of water delivery to the bed

Estimating Winter Balance and Its Uncertainty from Direct Measurements of Snow Depth and Density on Alpine Glaciers

Accurately estimating winter surface mass balance on glaciers is central to assessing glacier health and predicting glacier run-off. However, measuring and modelling snow distribution is inherently difficult in mountainous terrain. Here, we explore rigorous statistical methods of estimating winter balance and its uncertainty from multiscale measurements of snow depth and density. In May 2016, we collected over 9000 manual measurements of snow depth across three glaciers in the St. Elias Mountains, Yukon, Canada. Linear regression, combined with cross-validation and Bayesian model averaging, as well as ordinary kriging are used to interpolate point-scale values to glacier-wide estimates of winter balance. Elevation and a wind-redistribution parameter exhibit the highest correlations with winter balance, but the relationship varies considerably between glaciers. A Monte Carlo analysis reveals that the interpolation itself introduces more uncertainty than the assignment of snow density or the representation of grid-scale variability. For our study glaciers, the winter balance uncertainty from all assessed sources ranges from 0.03 to 0.15 m w.e. (5–39%). Despite the challenges associated with estimating winter balance, our results are consistent with a regional-scale winter-balance gradient

Pursuit of Optimal Design for Winter-Balance Surveys of Valley-Glacier Ablation Areas

Efficient collection of snow depth and density data is important in field surveys used to estimate the winter surface mass balance of glaciers. Simultaneously extensive, high resolution, and accurate snow-depth measurements can be difficult to obtain, so optimisation of measurement configuration and spacing is valuable in any survey design. Using in-situ data from the ablation areas of three glaciers in the St. Elias Mountains of Yukon, Canada, we consider six possible survey designs for snow-depth sampling and N = 6–200+ sampling locations per glacier. For each design and number of sampling locations, we use a linear regression on topographic parameters to estimate winter balance at unsampled locations and compare these estimates with known values. Average errors decrease sharply with increasing sample size up to N ≈ 10–15, but reliable error reduction for any given sampling scheme requires significantly higher N. Lower errors are often, but not always, associated with sampling schemes that employ quasi-regular spacing. With both real- and synthetic data, the common centreline survey produces the poorest results overall. The optimal design often requires sampling near the glacier margin, even at low N. The unconventional “hourglass” design performed best of all designs tested when evaluated against known values of winter balance

Detection and characterization of discontinuous motion on Thompson Glacier, Canadian High Arctic, using synthetic aperture radar speckle tracking and ice-flow modeling

We investigate unusual discontinuous glacier motion on Thompson Glacier, Umingmat Nunaat, Arctic Canada, using synthetic aperture radar (SAR) images and ice-flow modeling. A novel intensity-rescaling scheme is developed to reduce errors in high-resolution speckle tracking, resulting in a ~25% improvement in accuracy. Interferometric SAR (InSAR) and speckle tracking using high resolution RADARSAT-2 data indicate velocity discontinuities of up to 1 cm d−1 across deep and longitudinally extensive supraglacial channels on Thompson Glacier. We use a cross-sectional finite-element ice-flow model to determine the conditions under which velocity discontinuities of the observed magnitude and signature are possible. The modeling suggests that discontinuous motion across (long and straight) supraglacial channels can occur without ice fracture and under a wide variety of glacier thermal structures, including in fully temperate glaciers. Despite the wide range of conditions conducive to discontinuous motion, the form we observe requires that the associated channels be deep, longitudinally extensive and located in regions of lateral shearing. We speculate that these combined conditions are rare except on polythermal glaciers, where drainage features such as moulins are comparatively scarce and lower deformation rates allow channels to incise consistently and persist over many years

A stationary impulse-radar system for autonomous deployment in cold and temperate environments

L. Mingo acknowledges the SR&ED Tax Incentive Program of the CRA for offsetting some R&D costs of the sIPR. G. Flowers and D. Bigelow are grateful to NSERC, CFI, CSA, NSTP, PCSP and SFU for funding, Kluane First Nation, Yukon Government and Parks Canada for access to Yukon field sites and S. Williams, L. Goodwin, A. Pulwicki, J. Crompton,F. Beaud and Canadian astronaut D. Saint-Jacques for support and field assistance. C. Schoof and C. Rada provided time-lapse imagery and 2017 water-pressure data. A. Crawford and D. Mueller acknowledge funding from NSERC, Transport Canada, PKC/NSTP and the Garfield Weston Foundation, as well as support from Arctic Net, the CCGS Amundsen crew and pilots O. Talbot, A. Roy and G Carpentier.Stationary ice-penetrating radar (sIPR) systems can be used to monitor temporal changes in electromagnetically sensitive properties of glaciers and ice sheets. We describe a system intended for autonomous operation in remote glacial environments, and document its performance during deployments in cold and temperate settings. The design is patterned after an existing impulse radar system, with the addition of a fibre-optic link and timing module to control transmitter pulses, a micro-UPS (uninterruptable power supply) to prevent uncontrolled system shutdown and a customized satellite telemetry scheme. Various implementations of the sIPR were deployed on the Kaskawulsh Glacier near an ice-marginal lake in Yukon, Canada, for 44–77 days in summers 2014, 2015 and 2017. Pronounced perturbations to englacial radiostratigraphy were observed commensurate with lake filling and drainage, and are interpreted as changes in englacial water storage. Another sIPR was deployed in 2015–2016 on ice island PII-A-1-f, which originated from the Petermann Glacier in northwest Greenland. This system operated autonomously for almost a year during which changes in thickness of the ice column were clearly detected.Publisher PDFPeer reviewe

An Imbalancing Act: The Delayed Dynamic Response of the Kaskawulsh Glacier to Sustained Mass Loss

The Kaskawulsh Glacier is an iconic outlet draining the icefields of the St. Elias Mountains in Yukon, Canada. We determine and attempt to interpret its catchment-wide mass budget since 2007. Using SPOT5/6/7 data we estimate a 2007–18 geodetic balance of −0.46 ± 0.17 m w.e. a−1. We then compute balance fluxes and observed ice fluxes at nine flux gates to examine the discrepancy between the climatic mass balance and internal mass redistribution by glacier flow. Balance fluxes are computed using a fully distributed mass-balance model driven by downscaled and bias-corrected climate-reanalysis data. Observed fluxes are calculated using NASA ITS_LIVE surface velocities and glacier cross-sectional areas derived from ice-penetrating radar data. We find the glacier is still in the early stages of dynamic adjustment to its mass imbalance. We estimate a committed terminus retreat of ~23 km under the 2007–18 climate and a lower bound of 46 km3 of committed ice loss, equivalent to ~15% of the total glacier volume

Seismologically Observed Spatiotemporal Drainage Activity at Moulins

Hydrology is important for glacier dynamics, but it is difficult to monitor the subsurface drainage systems of glaciers by direct observations. Since meltwater drainage generates seismic signals, passive seismic analysis has the potential to be used to monitor these processes. To study continuous seismic radiation from the drainage, we analyze geophone data from six stations deployed at the Kaskawulsh Glacier in Yukon, Canada, during the summer of 2014 using ambient noise cross-correlation techniques. We locate the noise sources by backprojecting the amplitude of the cross correlation to the glacier surface. Most of the ambient noise sequences are found in two clusters, with each cluster located in the vicinity of a moulin identified at the surface. Stronger seismic radiation is observed during the day, consistent with expected variability in melt rates. We demonstrate that the sparse seismic network array with 2 km station separation has the ability to detect moulins within the array with a precision of 50 m. We confirm that seismic activity is correlated with air temperature, and thus, melt, on a diurnal timescale, and precipitation correlates with the activity at longer timescales. Our results highlight the potential of passive seismic observations for monitoring water flow into subglacial channels through moulins with an affordable number of seismic stations, but quantification of water flow rates still remains a challenge. The cross-correlation backprojection technique described here can also potentially be applied to any localized source of ambient noise such as ocean noise, tectonic tremor, and volcanic tremor

Seismologically Observed Spatiotemporal Drainage Activity at Moulins

Hydrology is important for glacier dynamics, but it is difficult to monitor the subsurface drainage systems of glaciers by direct observations. Since meltwater drainage generates seismic signals, passive seismic analysis has the potential to be used to monitor these processes. To study continuous seismic radiation from the drainage, we analyze geophone data from six stations deployed at the Kaskawulsh Glacier in Yukon, Canada, during the summer of 2014 using ambient noise cross-correlation techniques. We locate the noise sources by backprojecting the amplitude of the cross correlation to the glacier surface. Most of the ambient noise sequences are found in two clusters, with each cluster located in the vicinity of a moulin identified at the surface. Stronger seismic radiation is observed during the day, consistent with expected variability in melt rates. We demonstrate that the sparse seismic network array with 2 km station separation has the ability to detect moulins within the array with a precision of 50 m. We confirm that seismic activity is correlated with air temperature, and thus, melt, on a diurnal timescale, and precipitation correlates with the activity at longer timescales. Our results highlight the potential of passive seismic observations for monitoring water flow into subglacial channels through moulins with an affordable number of seismic stations, but quantification of water flow rates still remains a challenge. The cross-correlation backprojection technique described here can also potentially be applied to any localized source of ambient noise such as ocean noise, tectonic tremor, and volcanic tremor

Five decades of radioglaciology

Radar sounding is a powerful geophysical approach for characterizing the subsurface conditions of terrestrial and planetary ice masses at local to global scales. As a result, a wide array of orbital, airborne, ground-based, and in situ instruments, platforms and data analysis approaches for radioglaciology have been developed, applied or proposed. Terrestrially, airborne radar sounding has been used in glaciology to observe ice thickness, basal topography and englacial layers for five decades. More recently, radar sounding data have also been exploited to estimate the extent and configuration of subglacial water, the geometry of subglacial bedforms and the subglacial and englacial thermal states of ice sheets. Planetary radar sounders have observed, or are planned to observe, the subsurfaces and near-surfaces of Mars, Earth's Moon, comets and the icy moons of Jupiter. In this review paper, and the thematic issue of the Annals of Glaciology on ‘Five decades of radioglaciology’ to which it belongs, we present recent advances in the fields of radar systems, missions, signal processing, data analysis, modeling and scientific interpretation. Our review presents progress in these fields since the last radio-glaciological Annals of Glaciology issue of 2014, the context of their history and future prospects